Stability of Quantum Motion: Beyond Fermi-golden-rule and Lyapunov decay 
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We study, analytically and numerically, the stability of quantum motion for a classically chaotic 
system. We show the existence of different regimes of fidelity decay. In particular, when the 
underlying classical dynamics is weakly chaotic, deviations from Fermi- Golden-rule and Lyapounov 
regimes are observed and discussed. 
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The nature of correlations decay is an important sub- 
ject in different fields of physics. In particular, after the 
discovery of the so-called dynamical chaos, a large effort 
has been devoted to understand their behavior in relation 
to dynamical properties. The main reason is to know the 
precise conditions under which a statistical description is 
legitimate and to estimate the nature of the approxima- 
tions which are involved. 

Another important characteristic of dynamical systems 
is the stability of their solutions under slight variation of 
the Hamiltonian. A quantitative measure of this stability 
is given by the so-called fidelity or quantum Loschmidt 
echo. The fidelity M(t) — \m(t)\ 2 , measures the over- 
lap of two states started from the same initial state 
and evolved under slightly different Hamiltonians H and 
H = H n + eV, which are classically chaotic, 




m(t) = (^ \exp(iHt/h)e^p(-iH Q t/h)\^a} 



(1) 



Quite surprisingly, in spite of its physical relevance, 
the behavior of fidelity has been scarcely considered and 
only recently, in connection with quantum computation, 
a large number of papers appeared. Some important fea- 
tures of fidelity are now understood even though we are 
still far from the detailed level of knowledge we have 
about related quantities such as correlations functions 
and escape probabilities. So far, above the perturbative 
regime of small e, with Gaussian type decay Q, 0, 0, 
two main types of exponential decay of the fidelity have 
been identified: i)- The Fermi-golden-rule (FGR) decay, 
with the exponent given by the half-width of the cor- 
responding local spectral density of states^, 0, IE 01 • 
This decay has been related to the decay of autocor- 
relation function 0. ii)- The Lyapunov regime, above 
the FGR regime, with decay rate given by the Lya- 
unov exponent of the underlying classical dynamics 

"HIlllliaOEHElEi. 

In this paper we show that for classically chaotic sys- 
tems, in particular those with weak chaos, the behav- 
ior of fidelity can be much more rich and complex than 



FIG. 1: Fidelity M(t) as a function of a 2 t for K = 0.4, 
e « 7.67 x 10" 5 and N = No, 2N , 4iV , 8N (from bottom 
to top) where iV = 4096. (<r = 0.05,0.1,0.2,0.4). The FGR 



decay : 



is shown by the dashed line. Full circles rep- 



resent the semiclassical values M a (i) at a = 0.4, computed 
with expression (J7J. The numerically computed semiclassi- 
cal values Mf(t) turn out to be negligible so that M(t) is 
well approximated by M a (t), as clearly seen from the fig- 
ure at o = 0.4. Averages were performed over 400 initial 
point sources, with 9o taken randomly in the interval [0, 2n). 
(The same decaying behaviors are observed for initial Gaus- 
sian wavepackets.) 



expected. In particular we study perturbation borders 
which separate different types of decay. 

We start by displaying numerical results which 
strongly deviate from the expected behavior. We 
consider here the simple, well known, sawtooth map 
model^l- The classical map writes: 

p = p + K (8-Tr),6 = 6+p. (mod 2tt) (2) 

For Kq > 0, the motion is completely chaotic, with Lya- 
punov exponent A = ln{(2 + K + [{2 + K ) 2 -4] 1 / 2 )/2}. 
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FIG. 2: The exponential decay rate 7 versus perturbation 
strength a, calculated from the best fit of lnM(i). Gaussian 
wavepackets are taken as initial states. The solid curve shows 
the rate V ~ 2.2cr 2 of the FGR-decay. The dashed horizonal 
lines correspond to the Lyapunov exponents A = 0.62, 0.96 
and 1.32 for K = 0.4, 1 and 2, respectively. N = 131072. 



The quantum evolution on one map iteration is described 
by 

^=U ip = exp [-if/(2h)] exp[ik (6 - n) 2 /2}^, (3) 

where p = —ihd/dO and fco = Kq/H, with the effective 
Planck constant H = 2n/N and N being the dimension of 
the Hilbert space. For the perturbed system, K = K + e 
and k — fco + a, where a — e/H and e << Kq. 

In Fig. H we show the fidelity decay in the expected 
FGR regime 1/y/N < a < 1. In spite of the fact that 
the classical motion is chaotic, it is clearly seen that the 
behavior does not obey the FGR which, according to 
[1G3, should be cx exp{-Tt) with T » 2.2a 2 . The same 
conclusion can be drawn from Fig. where we plot the 
decay rate 7 of the fidelity as a function of a. Indeed at 
Kq = 0.4 the decay rate 7 versus a appears quite different 
from the quadratic one 16]. 

Deviations are present even at Kq = 1 and only at 
Kq = 2 one has good FGR decay. Moreover, above the 
FGR regime, where one expects Lyapunov decay, there 
are strong oscillations above and below the decay rate A 
(for Kq = 1, and 2). Only at larger a values, one enters 
the Lyapunov regime. 

In order to explain the above numerical results, we 
start from the standard semiclassical approach 0, 0|- 
For simplicity, we consider a finite configuration space, 
with dimension d and volume V — J dr. The momen- 
tum space is also finite, with a volume V p . In the semi- 
classical approach, an initial state ^>o( r o) is propagated 
by the semiclassical Van Vleck-Gutz wilier propagator, 
ijj(r;t) = J dr if S c(r,r ;i)V>o(ro), where K sc (r, r ; t) = 



K s (r,r ;t) 



1/2 



(2nih) d /' 



: exp 



% Z7T 

-5 s (r, r ;t) - 



(4) 



The label s in Eq. {3} (more exactly s(r, ro;t)), indi- 
cate classical trajectories starting at ro and ending at r 
in a time t; S s (r,ro;t) is the time integral of the La- 
grangian along the trajectory s, S s (r,r ;t) — J dt'C, 
C s = |dct(9 2 S s /droidrj)\, and /j, s is the Maslov index 
counting the conjugate points. 

In Ref. 0|> ^ is shown that the semiclassical approx- 
imation to m(t) for initial Gaussian wavepackets has a 
simple and convenient expression in the initial momen- 
tum space. Following similar arguments for initial point 
sources, (r|Wo) = \J (2nh) d /V p 8(r — r ), ( the theory can 
be extended to general initial states), one can write m(t) 
as 



m(r ,i) 



dpoexp 



: AS(po,r ;i) 



(5) 



where AS^po, i"o; £) is the action difference along the 
trajectory starting at (ro,po) for the two systems Hq 
and H. In the first order classical perturbation theory, 
AS , (po,r ;t) = e / dt'V[r(t')], with V evaluated along 
the trajectory. 

The averaged (over r ) fidelity can be separated into 
a mean- value part and a fluctuating part [llj, denoted 
by M a (t)_and Mf(t) respectively, M(t) = \m(t)\ 2 = 
M a (t) +Mf(t), where 



M a {t) ee |7fl(<)| 2 , with fn(t) 



dr m{ro,t). (6) 



From Eqs. an d © 5 it is seen that the mean- value 
part M a (t) can be expressed in terms of the distribution 
P(AS) of the action difference AS, 



M a (t) 



dASe lAS/h P(AS) 



where 



(7) 



P(AS) = , 1 - I drodpoS [AS - AS( Po ,r ; t)} . (8) 
/ dr dp J 



I 

It is usually assumed that for chaotic systems P(AS) 
is close to a Gaussian with a variance [2e 2 K(E)t], where 
K(E) = J a °° dt{V\r(t)]V[(r(0)]) is the classical action dif- 
fusion constant ,3j- As a result, M a (t) ~ e~ rt , where 
r = 2a 2 K{E). At small a, the fluctuation is small com- 
pared with the average value, because the phase on the 
right hand side of Eq. (|SJ) is proportional to a; then, 
M(t) ~ M a (t) has the FGR-decay. 

Let us now consider a fixed ro, and divide the space 
of the initial momenta po into connected, disjoint sub- 
spaces, denoted by A a , where each A a is the largest pos- 
sible subspace such that the correspondence between po 
and the final position r is one-to-one, i.e. different po 
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FIG. 3: Comparison between the exact M(t), its semiclassical 
mean- value part M a , calculated by using Eqs. © and ©, 
and the fluctuation part Mf = M sc — M a , where M BC is the 
semiclassical approximation to the fidelity, computed by the 
expression 0. Here N = 131072, K = 2^_and a = 0.9 
(left panel), a = 3(right panel). The exact M(t) is in good 
agreement with its semiclassical approximation M sc (t). The 
average is taken over 500 initial point sources. 



FIG. 4: _Similar to Fig. [2 for K = 1 anda = 6. At this 
large a, M a (t) is negligible compared with Mf(t). 
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inside each single component A a gives different final po- 
sition r. It is always possible to make such a division. 
The number of subspaces A a is denoted by N a . Note 
also that the sizes of A a decrease exponentially with in- 
creasing time t. When po runs over a subspace A a , r 
may run over part of the configuration space, denoted by 
V a . Note that, with this division of the po subspace, the 
trajectories starting at ro are divided into N a groups and 
"near" trajectories typically belong to the same group. 

The amplitude m(ro,i) in Eq. J5j can now be written 
as m(r ,t) ~ J] Q TO a( r o,*), where 



m a (r ,t) 



1 



drC s exp 



-AS s (r,r ;i) 



(9) 



with integration over the subspace V Q , in which the 
change of variable po — » r within the subspaces A a has 
been done and AS s (r, r ;i) coincides with AS^po, r ; t) 
for the same trajectory s starting at (ro,po) with po € 
A a - M f(t) is written as 



M f {t) 



with m a f = m a (ro 1 1) ~ 



m(t) 



(10) 



When a is large enough, above a critical border <r/, 
m a (vo,t) can be regarded as possessing random phase, 
and therefore Mt can be approximated by its diagonal 



dr oE 



drC s exp ( xASf 



,(H) 



where the second approximation is obtained by noticing 
that \m(t)/N a \ << \m a \ at large a. 

When the phase space is homogeneous with constant 
local (maximum) Lyapunov exponent A, as in the saw- 
tooth map, the number N s (ro, r) of trajectories connect- 
ing two points ro and r in the configuration space is ap- 
proximately N s ~ N a ~ e xt |l7j. The summation over a 
in llljl gives a contribution approximately proportional 
to N s . At t large enough, the main time dependence 
of | J v drC s e^ ASs \ is given by C s cx e~ xt . Combining 
these results, it is seen that at a > <jf, Mf(t) has the 
Lyapunov decay, Mf(t) oc e _At . 

In order to have the Lyapunov decay for M(t), the 
term M a (t) must be small. To this end one needs to 
further increase a above a critical value ay, so that the 
variance of the phase of m(ro,t) with respect to ro will 
become so large that M a (t) is negligible. 

The right panel of Fig. Ogives an example of M a (t) w 
M f (t) . This explains the fluctuation of 7 versus a shown 
in Fig. |21 at K — 2 and a < 3. Fig. ^instead gives an 
example with a large enough (a > oy), so that M a (t) is 
negligible and M(t) ~M f (t). 

The deviation from FGR decay observed in Figs.n] an d 
El is due to the deviation of P(AS) from the Gaussian 
behavior. Indeed, when chaos in the underlying classical 
dynamics is strong enough (Kq > 1), correlations be- 
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FIG. 5: Distribution P[(AS- (AS))/e] of the classical action 
difference AS, at t — 10, calculated by taking randomly 10 7 
initial points in the phase space, where (AS) = et{V(6)) = 
— 7r 2 ei/6, with an average over the phase space. 



tween non-overlapping parts of a trajectory decay very 
rapidly and the distribution P(AS/e) reaches, in a rela- 
tively short time, the Gaussian distribution. This is the 
case of Fig.Elfor K = 2 , where K(E) = vr 4 /90 ~ 1.08 
with T = 2K(E)a 2 ~ 2.16a 2 , in agreement with the nu- 
merical results in Ref. |l2( and in Fig. [21 However, when 
Kq is not sufficiently large, e.g., Ko = 0.4, a consider- 
able deviation of P(AS/e) from the Gaussian distribu- 
tion appears for times comparable to the fidelity decay 
times (Fig. EJ|. According to Eq. 10, this leads to de- 
viations from the FGR-decay as observed in Fig. ^ Wc 
would like to draw the reader attention to the fact that, 
for Kq < 1 the saw-tooth map, even though completely 
chaotic, it possess a structure of cantori which, in the 
quantum case, can act as perfect barriers to quantum 
motion thus leading to localization of wavefunctions. 

Notice that the deviation of P(AS/e) from the Gaus- 
sian distribution depends on Kq but not on e or a. There- 
fore, by increasing a, the effect of this deviation becomes 
more and more important, since the FGR exponential 
decay has a decay rate proportional to a 2 while the devi- 
ation from the Gaussian remains unchanged. Therefore, 
for a given system, there is a critical value aa, below 
which the FGR decay is obeyed with good accuracy and 
above which FGR breaks down. This case is illustrated in 
Fig.[21for the case Ko = 1, which coincides with the well- 
known Arnold cat map, the paradigmatic model of chaos. 
Here the distribution P(AS/e) (at t= 10) is slightly dif- 
ferent from the Gaussian distribution and the decay rate 
7 of fidelity deviates from the FGR decay for a > 0.3. 
In cases of weak classical chaos, the value of ad can be 
so small that FGR is never observed (e.g. the case with 
Ko = 0.4). The left panel in Fig. [3]shows instead a case 
at K o = 2 and a = 0.9, in which M a (t) obeys the FGR 



decay and Mf(t) is negligible. 

To summarize: above the perturbative border, the fi- 
delity has a FGR decay for a < ad-, while for a > a r , 
it has the Lyapunov decay. In the intermediate region, 
for ad < a < af, the fidelity deviates from FGR and 
can decay even faster than Lyapunov. For <jf < a < a r , 
M a {t) ~ Mf{t) and the decay rate of M{t) fluctuates 
around the Lyapunov exponent. It may be useful to recall 
here the physical meaning of different borders. Above ad, 
the distribution P(AS) deviates from the Gaussian and 
this induces deviations from the expected FGR decay. 
Below a r , M a (t) is non negligible as compared to M*{£) 
and this induces deviations from the expected Lyapunov 
decay. 

It may be interesting to remark that the relation be- 
tween the decomposition in two part of M(t) here and 
that in Ref. Q is the following. At a small enough 
M(t) ~ M a {t) ^ M nd (t), with_M/(t) and M d {t) neg- 
ligible; while at a large enough M(t) ~ M f (t) ~ M d {t), 
with M a {t) and M nd (t) negligible. In the intermediate 
regime of a, in particular, in the crossover from the FGR 
decay to the Lyapunov decay, there may be considerable 
difference between the two divisions. 

In this paper, by using the sawtooth map, we have 
demonstrated that the fidelity decay in a generic chaotic 
system can have a very complex behavior. In partic- 
ular, deviations from the Fermi golden rule (for weak 
chaos) and Lyapunov decay have been discussed as well 
as the existence of perturbation borders separating differ- 
ent regimes. It is our opinion that fidelity is an important 
quantity which characterizes the stability of classical and 
quantum systems. It therefore deserves deeper analyti- 
cal and numerical studies in order to fully understand its 
behavior in different dynamical regimes. 
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